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^ 1 Abstract: An analysis of the gravitational instability in presence of dissipa- 

(3JT)! tive effects is addressed. In particular, the standard Jeans Mechanism and the 

generalization in treating the Universe expansion are both analyzed when bulk 
\ viscosity affects the first-order Newtonian dynamics. As results, the perturbation 

■ evolution is founded to be damped by dissipative processes and the top-down 

mechanism of structure fragmentation is suppressed. In such a scheme, the value 
of the Jeans Mass remains unchanged also in presence of viscosity. 

O 

o 

O ■ 1 Introduction 
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The Cosmological Standard Model pQ well describes many parts of the Universe evolution 
and it takes into account the Friedmann-Lemaitre-Robertson- Walker (FLRW) metric as the 
highest symmetric background. In fact, considering the mean energy-density at big scales, 
i.e., greater than 100 Mpc, it tends to an homogeneous distribution. On the other hand, the 
Universe observation at small scales shows a very inhomogeneous and anisotropic matter- and 
energy-distribution. In this respect, the isotropic hypothesis of the Cosmological Principle 
2\ is not based on the big-scale observations but on the strong isotropy of the Cosmic 
Background Radiation, which has a black-body spectral-distribution at T ~ 2.73 K with 
temperature fluctuations of order O(10 -4 ). Moreover, the Cosmological Standard Model is 
confirmed by the primordial- nucleosynthesis prediction for the light elements, which is in 
agreement with direct observations. 

Such a dichotomy between the isotropy of region at red-shift z rs ~ 10 3 and the extreme 
irregularity of the recent Universe, z rs <C 1, is at the ground of the interest in the study 
of the gravitational instability for the structure formation. The study of the cosmological 
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perturbation evolution can be separated in two distinct regimes: the linear regime, concerning 
density contrast 5 much less than the unity, and the non- linear one in which <5 > 1, giving 
rise to the effective structure formation. Despite the approximate hypotheses, the linear 
regime provides interesting predictive information also at low red-shift, since an analytical 
description can be addressed to study the growth of the density contrast. 

As matter of fact, we underline that the study of the perturbation dynamics in the 
radiation-dominated Early Universe requires a pure relativistic treatment, in order to cor- 
relate the matter fluctuations with the geometrical ones H|. On the other hand, the 
evolution during the matter-dominated era can be consistently described using the Newto- 
nian approximation picture, as soon as sub-horizon-sized scales are treated. In this scheme, 
the fundamental result of the density-perturbation analysis is the so-called Jeans Mass, which 
is the threshold value for the fluctuation masses to condense generating a real structure. If 
masses greater than the Jeans Mass are addressed, density perturbations begin to diverge as 
function of time giving rise to the gravitational collapse [3 [S] . 

In this work, we are aimed to consider dissipative effects into the fluid dynamics, in the 
linear Newtonian regime. The starting point is the Eulerian set of equations governing the 
fluid motion (Section 2) on which one can develop a perturbative theory by adding small 
fluctuations to the unperturbed background (see also pJ[B])- In particular, we introduce in 
the first-order analysis the so-called bulk viscosity (we neglect the shear or first viscosity since 
we are dealing with homogeneous model and no internal frictions arise). Such kind of viscosity 
can be expressed in terms of the thermodynamical parameters of the fluid. Following the line 
of the "Landau School" [9] [TO] [HJ [12l [13] , we phenomenologically describe this quantity as 
a function of the Universe energy density p via a power-law: £ = z p s where s = const, and 
z is a parameter which defines the intensity of viscous effects, see also [T4l [T51 fTB] . 

Two different cases are treated: the standard Jeans Mechanism [5J [2] and the generaliza- 
tion in presence of the expanding Universe [T7j. In the first approach (Section 3), the 
unperturbed background is assumed to be characterized by a static and uniform solution 
of the fluid parameters while, in the second case (Section 4), we consider the effects of the 
Universe expansion on the zeroth-order dynamics which is now described by the Friedmann 
Eq. 's for an homogeneous and isotropic Universe. 

As a result, we show how the presence of bulk viscosity damps the density-contrast growth, 
suppressing the structure formation in the top-down fragmentation scheme pQ (mainly as- 
sociated with the hot dark-matter phenomenology [HI [HI [20l EI]) without changing the 
threshold value of the Jeans Mass. In fact, in both cases, the grater the viscosity coefficient 
is, the slower the density contrast diverges in time, in the limit t — > oo. In this respect, con- 
sidering a collapsing macro-structure, i.e., of mass greater than the Jeans Mass, if viscosity 
becomes sufficiently large then the Jeans fragmentation mechanism is deeply unfavored since 
the sub-scale density-contrasts are damped. In particular, in treating the Jeans Model, bulk 
viscosity affects the perturbation evolution generating a new decreasing regime in place of 
the pure oscillatory behavior addressed in the Jeans analysis. 

The main merit of this work is to be determined in having traced a possible scenario 
for fragmentation processes in presence of viscosity. We infer that the unfavored nature of 
the top-down mechanism, appearing when a viscous trace is present, can survive also in 
the non-linear regime when dissipative effects play surely an important role in the structure 
formation. 

2 Motion Equations of Viscous Fluids 

In order to describe the Newtonian evolution of a fluid, we here want to introduce the Eulerian 
equations governing the fluid parameters: the density p, the local 3- velocity v (of components 
v a ) and the pressure p, in presence of a gravitational potential <j>. 
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Adiabatic ideal fluids are governed, in Newtonian regime, by the following set of equations 
[22j : the Continuity Eq., which guarantees the energy conservation 

d t p + V-(pv) = 0, (1) 

the Euler Eq., which ensures the momentum conservation 

pftv + p(v-V)v = -Vp-pV0, (2) 

while pressure and density are linked by the Eq. of State (EoS): p = p (p). In this picture, 
the sound speed is defined by the relation v 2 s = 6p/Sp. 

Let us introduce the effects of the energy dissipation during the motion of the fluid, due 
to the thermodynamical non-reversibility and to internal friction (we neglect the thermal 
conductivity). To obtain the motion equations for a viscous fluid, we have to include some 
additional terms in the ideal fluid description. The Continuity Eq. is derived by the time 
evolution of the matter density and by the mass conservation law. This way, it remains 
valid for any kind of fluid. Euler Eq., in absence of the gravitational field, rewrites (here 
a = 1,2,3) 

d t {pv a ) = -dpllap , (3) 

where n Q( g denotes the momentum-flux energy-tensor. If ideal fluids are addressed, we deal 
with completely reversible transfer of momentum, obtaining the expression: II Q/ 3 = p5 a j3 + 
pv a vp. Viscosity is responsible for an additional term d a p to this expression due to another 
irreversible momentum transfer, where non-vanishing velocity gradients are present. For a 
viscous fluid we get 

n a/ 3 =pS a f3 + pv a vp - a a p = -a af 3 + pv a V/3 , a al 3 = -p5 a p + a af 3 , (4) 

where a a p is the stress tensor and a a p is called the viscous stress-tensor. 

The general form of a a p can be derived by a qualitative analysis of the velocity gradients 
in presence of uniform rotation and volume changes of the fluid. The most general form of 
the viscous stress tensor is [22] 

cr«/9 = f] (dpv a + d a V/3 - | <5 Q/ 3<9 7 u 7 ) + C S a p d~fV~f , (5) 

where the coefficients ?j e ( are not dependent of velocity (the fluid is isotropic and its prop- 
erties must be described only by scalar quantities) and the term proportional to rj coefficient 
vanishes for the a and (3 contraction. Here, the coefficient r\ is called shear viscosity while £ 
denotes bulk viscosity and they are both positive quantities. 
Using Continuity Eq., the ideal fluid Euler Eq. rewrites 

p(d t v a + dpVa) = -d a p , 

and the motion equation of a viscous fluids can now be obtained by adding the expression 
d[3<T a p to the rhs of the equation above, obtaining 

p(d t v a + vp dpv a ) = -d a p + dp[r) (dpv a + d a vp - § <5 Q/3 (9 7 i; 7 )] + <9 Q (C <9 7 -u 7 ) . 

The viscous coefficients are not constant and we have to express their dependence on the 
state parameters of the fluid. Since we are interested to treat isotropic and homogeneous 
perturbative cosmological models, wc can safely neglect the first viscosity (shear viscosity) in 
the unperturbed dynamics. In fact, in such models there is no displacement of matter layers 
with respect to each other and this kind of viscosity represents the energy dissipation due 
to this effect. Indeed, in presence of small inhomogeneities such affect should be taken into 
account, in principle. However, in this work, we are aimed at studying the behavior of scalar 
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density perturbations. In this respect, volume changes of a given mass scale are essentially 
involved and, therefore, we concentrate our attention to bulk-viscosity effects only. In fact, 
we expect that the non-equilibrium dynamics of matter compression and rarefaction is more 
relevant than friction among the different layers. 

According to literature developments UHl HU HH HE] , we now assume the bulk- viscosity 
coefficient as a function of the energy density p(t) expressed via a power-law of the form 

( = zp s , (6) 

where s = const, and z is a constant parameter which defines the intensity of viscous effects. 
With these assumptions, the Euler Eq. takes the following form 

p9 t v + p(v- V)v + Vp -CV(V-v) = 0, (7) 

which is the well-known Navier-Stokes Eq. 

This analysis is developed without considering the gravitational field, which has to be 
introduced in the Euler Eq. as usual. We have also to consider the equation describing the 
gravitational field itself: the Poisson Eq. Let us now recall the set of motion equations in 
the case of an adiabatic viscous fluid: 

d t p + V • (pv) = , (8a) 
p d t v + p (v • V)v + Vp - C V(V ■ v) + p V0 = , (8b) 

V 2 - AirGp = , (8c) 

such a system is the starting point to analyze the gravitational instability. 

3 Analysis of the dissipative Jeans Mechanism 

The Universe is uniform at big scales but many concentrations are presented at small scales, 
e.g., galaxies and clusters, where the mass density is larger than the Universe mean-density. 
These mass agglomerates are due to the gravitational instability: if density perturbations are 
generated in a certain volume, the gravitational forces act contracting this volume, allowing 
a gravitational collapse. The only forces which contrast such gravitational contraction are 
the pressure ones, which act in order to maintain uniform the energy density. The Jeans 
Mechanism analyzes what are the conditions for which density perturbations become unstable 
to the gravitational collapse. 

This model [5] is based on a Newtonian approach and the effects of the expanding Universe 
are neglected. The fundamental hypothesis of such an analysis is a static and uniform solution 
for the zeroth-order dynamics 

v = , Po — const. , p = const. , <p = const. (9) 

Of course, this assumption contradicts the gravitational equation, but we follow the original 
Jeans analysis imposing the so-called "Jeans swindle" We underline that our study will 

focus on Universe stages when the mean density is very small: in particular the recombination 
era, after decoupling. This way, the effects of bulk viscosity on the unperturbed dynamics 
can be consistently neglected in view of its phenomenological behavior (JHJ) . 

3.1 Remarks on the Jeans Model 

In the standard Jeans Model, a perfect fluid background is assumed. After setting ( = 
in the motion equations ||5J), let us add small fluctuations to the unperturbed solution: p — 
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p + Sp, p = p + Sp, 4> = (f) + S(j), v = v + 5v. Furthermore, only adiabatic perturbations are 
treated and the sound speed is defined as v 2 = dp/ dp. Substituting such expressions in the 
system ©, after standard manipulation (we remark that second-order terms are neglected), 
one differential equation for the density perturbations can be derived: 

d 2 Sp - v 2 V 2 Sp = AkG p Sp . (10) 

To study the properties of Sp, we now consider a plane-wave solutions of the form 

Sp(r,t) = Ae lut - l ^ r , (11) 

where w and k (fc = |k|) are the angular frequency and the wave number, respectively. This 
way, one can obtain the following dispersion relation 

lu 2 = v 2 s k 2 -AnGp a . (12) 

In this scheme, two different regimes are present: if u 2 > a pure time oscillatory-behavior 
for density perturbations is obtained. While if u 2 < 0, the fluctuations exponentially grow 
in time, in the t — * oo asymptotic limit (i.e., we choose the negative imaginary part of the 
angular frequency solution) and the gravitational collapse is addressed since also the density 
contrast S = Sp/p diverges. The condition uj 2 — defines the so-called Jeans Scale Kj and 
the Jeans Mass Mj (which is the total mass in a sphere of radius R = ir/Kj). Such threshold 
quantities read 

Let us now analyze in some details the two regimes. In the case M < Mj (i.e., to 2 > 0), 
Sp behave like two progressive sound waves (with constant amplitude) propagating in the 
±k directions with velocity v w — v s (l — (Kj/k) 2 ) 1 ^ 2 . In the limit k — > oo, the propagation 
velocity approaches the value v s , and fluctuations behave like pure sound waves. On the 
other hand, if fc — > Kj, stationary waves are addressed (i.e., v w = 0). 

In the case M > Mj (i.e., uj 2 < 0), density perturbations evolve like stationary waves 
with a time dependent amplitude. In particular, choosing the negative imaginary part of 
the solution for oj, the wave amplitude exponentially explodes, generating the gravitational 
collapse. 

3.2 Jeans Mechanism in presence of bulk viscosity 

Let us now analyze how viscosity can affect the gravitational-collapse dynamics. We recall 
that the only viscous process which can be addressed in an homogeneous and isotropic model 
is bulk viscosity. As discussed at the beginning of this Section, we are able to neglect such 
kind of viscosity in the unperturbed dynamics, which results to be described by the static 
and uniform solution 

We now start by adding the usual small fluctuations to such a solution, i.e., Sp, Sp, 5<j), 5v. 
In treating bulk-viscosity perturbations, we use the expansion C = Co + where 

Co = C(po) = zpt = cmist. , 5( = Sp (d(/dp) + ... = zsp^Sp + ... . (14) 

Substituting all fluctuations in the system (|8|) , we get the first-order motion-equations of the 
model 

d t Sp + p V • Sv = , (15a) 
Pa d t Sv + v 2 VSp + Po VS<p - Co V(V ■ <Sv) = , (15b) 

V 2 S(I)~AitGSp = . (15c) 



5 



With some little algebra, one can obtain an unique equation for density perturbations, de- 
scribing the dynamics of the gravitational collapse: 

p Q d 2 t 5 P - p v 2 V 2 5p - C„ V 2 d t Sp = AkGpI Sp . (16) 

Using the linearity of the equation above, a decomposition in Fourier expansion can be 
performed. This way, plane waves solutions can be addressed, obtaining a generalized 
dispersion relation 

p Q uj 2 - i Co k 2 uj + p (4nGp - v 2 k 2 ) = . (17) 

As in the standard Jean Model, the nature of the angular frequency is responsible of two 
different regimes for the density-perturbation evolution. The dispersion relation has the 
solution 

W = i|^±Vtf, Q = -^+v 2 k 2 -AnGp , (18) 
2 p 

thus we obtain the time exponential-regime for uj ^ and a damped oscillatory regime for 
Q > 0. It's worth noting that the pure oscillatory regime of the ideal fluid Jeans Mechanism 
is lost. The equation uj = admits the solutions K\ and K 2 which read 

Kl = ^lf lT Jl-(^) 2 )*, Kl ,K 2 >0, Kl <K 2 . (19) 

2 Co V V V PoVs > ' 

The existence of such solutions gives rise to a constraint on the viscosity coefficient: Co ^ 
C c = PoV s /Kj. An estimation in the recombination ergQ after decoupling, yields to the 
value Cc = 7.38 ■ 10 4 gem' 1 s' 1 and confronting this threshold with usual viscosity (e.g., 
^Hydr. _ g ^ _ ±Q-7g cm -i s _1 ), we can conclude that the range Co ^ Cc is the only of physical 
interest. Finally we obtain: uj ^ for k ^ K\, K 2 ^ k and Q > for K\ < k < K 2 . 

Let us now analyze the density-perturbation exponential-solutions in correspondence of 
uj «: 0: 

5p~e mt , » = -^TV=^. (20) 

2po 

To obtain the structure formation, the amplitude of such stationary waves must grow for 
increasing time. The exponential collapse for t — > 00 is addressed, choosing the (+) sign 
solution, only if w > 0, i.e., k < Kj with Kj < K-y < K 2 . As a result, we show how the 
structure formation occurs only if M > Mj, as in the standard Jeans Model. 

The viscous effects do not alter the threshold value of the Jeans Mass, but they change 
the perturbation evolution and the pure oscillatory behavior is lost in presence of dissipative 
effects. In particular, we get two distinct decreasing regimes: for K\ < k < K 2 (i.e., uj > 0), 
we obtain a damped oscillatory evolution of perturbations: 

Sp ~ e ~ 2 "a 1 cos(v / C t) , (21) 

while, for Kj < k < K\ and K 2 < k, density perturbations exponentially decrease as 
Sp ~ e wt , with w < 0, in the limit t — > 00. 



lr The parameters are set as follows: the usual barotropic relation p = c 2 / 'p 1 1 is assumed and the 
constant p can be derived from the expression expression Mj l|13jl . Universe is dominated by matter and 
we can impose the values: Mj ~ 1O 6 M , 7 = 5/3, p c = 1.879 h? ■ 1CT 29 gcm~ 3 , h = 0.7, z = 10 3 
and p a = p c z 3 = 0.92 ■ 10 — 20 gcm~ 3 . Using these quantities one finds p = 9.034 ■ 10 -7 gcm~ i , 
v a = 8.39 ■ 10 cms and the threshold value f c . 
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3.3 Implication for the top-down mechanism 

As shown above, since the pure oscillatory regime does not occurs, we deal with a decreasing 
exponential or a damped oscillatory evolution of perturbations. This allows to perform a 
qualitative analysis of the top-down fragmentation scheme pQ, i.e., the comparison between 
the evolution of two structures: one collapsing agglomerate with M ^> Mj and an internal 
non-collapsing sub-structure with M < Mj. If this picture is addressed, the sub-structure 
mass must be compared with a decreasing Jeans Mass since the latter is inversely proportional 
to the collapsing agglomerate background mass. This way, as soon as such a Jeans Mass 
reaches the sub-structure one, the latter begins to condense implying the fragmentation. 
In the standard Jeans Model, this mechanism is always allowed since the amplitude for 
perturbations characterized by M < Mj remains constant in time. On the other hand, the 
presence of decreasing fluctuations in the viscous model, requires a discussion on the effective 
damping and an the efficacy of the top-down mechanism. Of course, such an analysis contrasts 
the hypothesis of a constant background density, but it can be useful to estimate the strength 
of the dissipative effects. 

We now study two cases for different values of the bulk- viscosity coefficient: („ < 1 and 
Co > 1. In this analysis, a perturbative validity-limit has to be set: we suppose 5p/p a ~ 0.01 
as the limit of the model and we use the recombination era parameters °, in particular the 
initial time of the collapse is define as the beginning of the matter-dominated Universe, i.e., 
t a = t MD = 1.39- 10 13 s. 

In correspondence of a very small viscosity coefficient (Fig.l), we consider a decreasing 
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Figure 1: Case Co = gem,- 1 s _1 . 

structure density contrast 83 



Galaxy density contrast: 5g 
Ms = 10 Mq - (normal line). 



Mg = 1O 12 M - (dashed line). Sub- 



Structure of mass Ms — 10 Mq within a collapsing galaxy with mass Mq = 10 12 Af Q , the 
Jeans Mass is Mj = 10 6 M©. The sub-structure wave-number K$ is in the region K\ < K$ < 
K2 and the density perturbations evolve like (|2"Tj) . Fluctuations have to be imposed small 
at the initial time t , this way, we consider density contrasts (Sq for the galaxy and 5s for 
the sub-structure) of 0(1O~ 3 ). In this scheme, the galaxy starts to collapse and the validity 
limit is reached at t* — 6.25 • 10 13 . As a result, in Fig.l we can show how the sub-structure 
survives in the oscillatory regime during the background collapse until the threshold time 
value t* . Thus, we can conclude that, if the viscous damping is sufficiently small, the galaxy 
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formation occurs. 

Let us now discuss the case £ > 1 (Fig. 2) by changing the sub-structure mass, which is 
now Ms = Mq. Here, the viscosity coefficient is greater than one and the damping effects is 
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Figure 2: Case Co = 14 gem 1 s . Galaxy density contrast: Sg - Mq — 10 12 Mq - (dashed line). Sub-structure 
density contrast 5 s - Ms — Mq - (normal line). 



stronger, fn fact, when the galaxy density contrast reaches the threshold value 6q — 0.1, we 
obtain 5s = I CP 5 . The top-down mechanism for structure formation results to be unfavored 
by the presence of strong viscous effects: the damping becomes very strong and the sub- 
structure vanishes during the agglomerate evolution. 



4 Generalization to the expanding Universe background 

Let us now analyze the behavior of small perturbations, using Newtonian equations, on the 
expanding-Universe background 2 , 17, 23 . The equations which describe the homogeneous- 
and isotropic-Universe evolution are the well-known Friedmann Eq.'s. Such equations are 
derived by Einstein Eq.'s using a perfect fluid energy- momentum tensor as the matter source 
of the gravitational held. It's worth remarking that, we can safely address a Newtonian 
scheme for astrophysical models, as soon as we treat problems in which the energy density 
is dominated by non-relativistic particles and in which the linear scales involved are small 
compared with the characteristic scale of the Universe. Such a matter-dominated Universe 
is described by the FLRW metric 

ds 2 = dt 2 - a 2 (t)d£ 2 , (22) 

with an EoS so that p ~ (p C /)), here a(t) is the scale factor of the Universe. 

As in the Jeans Model, we here use the power-law ([6]) to describe bulk viscosity, which 
results to be proportional to a positive power of the energy density in the matter-dominated 
scheme we consider. This way, being the matter density very small, we can consistently 
neglect viscosity in the unperturbed dynamics. 

The zeroth-order solution corresponds to the evolution of an homogeneous and isotropic 
Universe filled with the following source: T^ u = diag[p, — p, —p, —p]. The dynamics equa- 
tions are the energy- momentum conservation-law T" v = (for \x = 0), written in a co-moving 
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frame (here and in the following, the dot ( ' ) denotes the total derivative wrt time), 



Pa 
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P + 3- (p + p)=0, (23) 
a 

and the cosmological equation 

« 2 +/C = ^pa 2 , (24) 

where K, = const, is the curvature factor. In this picture, the unperturbed solutions are, 
setting p = 0, 

'(-) 3 , v = r-, V0 o = tr7rGp o , (25) 

where p and a are dimensional constants, r (r = |r|) denotes the radial coordinate vector 
and, of course, a(t) satisfies ([24|) . The solutions v and </> are derived from the Continuity 
Eq. (|5a|) and the Poisson Eq. (f8"c| respectively, while the Navier-Stokes Eq. results to be 
satisfied since the Friedmann Eq.'s hold. 

To obtain the time dependence of the parameter involved in the model, we limit our 
analysis to early times, since fluctuations arise from the recombination era and, furthermore, 
the Jeans Mass is so small for recent times that it is of little interest [17]. This way, the 
study is restricted to scale factors satisfy the condition a(t) <C a , so that a 2 , 8irpa 2 /3 S> 1 
and we can use the zero-curvature solution without loss of generality. Setting /C = in the 
cosmological equation ([2"4"|l and using the solution p (j2l)|) , one can get the following time 
dependence 

a~t*'\ Po = ^. (26) 

The study of the gravitational instability is characterized by the evolution of the density 
contrast and, in particular, of the small fluctuations. In this respect, we underline that 
v 2 = Sp/5p takes account for first-order terms and we have to explicitly write its time 
dependence during the Universe expansion. For a general specific heat ratio 7, we can 
assume that the pressure varies as pi and the speed of sound is find to be 

v s ~ t 1 -"* . (27) 

Such solutions characterize the background dynamics of the expanding Universe. It's 
worth noting that, in this generalization, the unperturbed dynamics is now a real solution of 
the zeroth-order equations and we do not have to apply the "Jeans swindle" static-solution 
assumption. 



4.1 Review of the non-dissipative case 

We want now to study the behavior of the density contrast, without the presence dissipative 
effects. Since we consider small scales, i.e., r <C a (r/a = 0), as the fluid motion-equations 
one can assume the system of Newtonian equations ([8]). In particular, we neglect the viscous 
term in equation (|8bp and we perform the usual perturbation theory. 

The resulting first-order motion-equations are spatially homogeneous [2] so one can address 
plane-wave solutions of the form 

Sp(r,t)=p 1 (t)e i ^ , 5v(r,t) = v 1 (t)e^ a , ty(r,t) = ^(i) . (28) 

The factor l/a(t) represents the wave-length reduction dues to the Universe expansion: q = 
|q| is the co-moving weave number, being k = q/a the physical one. To complete our analysis, 
in the limit r/a = 0, it is convenient to decompose the time depending velocity fluctuations 
Vi into two part: one transversal and one parallel to the q direction, respectively: 

Vl (t) =v^ + ^q e , q-v^=0, e = -4,(q. Vl ). (29) 
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It is also useful to express p x in terms of the density contrast: pi(t) = p S. 

A simple algebraic analysis of the first-order dynamics shows that two different types 
of normal modes arise. The Rotational Modes are described by vj 1 and simply decay as 
v^(i) ~ 1/a during the Universe expansion. On the other hand, the Compressional Modes 
are characterized by e e S and require a more interesting analysis. Such modes are described 
by the equation 

6 + — 6 + ( ^f- - 4ttG Po ) 5=0, (30) 

which reduces to the Jeans dispersion-relation (fT2"j) as soon as we set a = const, and consider 
the physical wave number k. Taking into account the zeroth-order time dependence (|26p and 
(f2"?|) . one finds that the density-contrast solution of equation (|3T)|) involves Bessel functions. 
Such special functions have different behavior corresponding to small or large proper argu- 
ment. If the argument is large, i.e., much greater than one, density contrast oscillate, on the 
other hand, it evolves like 

S ~ t-V6±5/6 ; ( 31 ) 

as soon as the Bessel argument is much less than unity. The condition which separates the 
two regimes, implying the gravitational collapse in the limit t — > oo (of course choosing the 
positive solution for 6), can be write as v 2 q 2 /a 2 < 6irGp , which is substantially the same as 
the Jeans condition derived from (fT2"| . It worth underling that the standard Jeans condition 
is perfectly recast if the parameter 7 of ([2"T|) is set to the value 4/3. 

In conclusion, we can infer that the dynamics of the expanding Universe does not modify 
(substantially) the value of the Jeans Mass which remains the threshold to address the 
gravitational collapse of structures. 



4.2 Bulk-viscosity effects on the density-contrast dynamics 

We are now aimed at introducing bulk-viscosity effects into the dynamics. As discussed above, 
such a dissipative effect can be consistently neglected from the zeroth-order analysis, since 
we are dealing with a matter-dominated Universe. This way, the unperturbed background 
on which develop the perturbative theory corresponds to the solution (125p of a Friedmann 
Universe. 

Adding small fluctuations to the Newtonian system ([8|) and neglecting second-order terms, 
we get the following set of equations 

d t Sp + 3-5p+-(r- V)5p + p V • <5v = , (32a) 
a a 

Po d t Sv + p -5v + po - (r • V)6v + v 2 s V Sp + p V 6<p ~ Co V (V • <5v) = , (32b) 
a a 

V 2 S(f> - AnGSp = , (32c) 

where the relation 8p — v 2 8p has been used and we recall that Co = C(aO> see (fPfj). As 
in the non-dissipative case, the plane-wave expansion (|28p for the fluid parameters can be 
addressed. Using the hypothesis r/a ~ 0, the system above reduces to: 

p 1 + 3-p 1 + ^(q-v 1 ) = 0, (33a) 
a a 

v 1 + -v 1 + ^-qp 1 -47riGap 1 -^ + 4 2 -q(q- v i) = °- ( 33b ) 
a ap q a z p„ 

Let us now follow the standard analysis and use the decomposition (129p in order to compare 
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our results wrt the non-dissipative ones. We finally get 



(- 



Pa a 2 



a j 
a 



= 



S-^-e =0 
a 



(34a) 
(34b) 
(34c) 



The Rotational Modes are not affected by viscosity. In fact, they are governed by equation 
(|34ap which has the solution 

v^(t) ~ 1/a , (35) 

as in the standard analysis presented above. On the other hand, the Compressional Modes 
are influenced by the presence of viscosity. In particular, combining together (|34bj) and (|34c|) . 
we get an equation which generalizes the compressional equation (|30[) . It reads 



/ a 
V a 



Col 

Po a* 



2 2 



AnGpo 5=0 



(36) 



This is the fundamental equation which governs the evolutions of the density contrast on 
an expanding Universe. Let us now write explicitly the time dependence of the parameters 
involved in the model. The zeroth-order analysis still remains valid in presence of viscosity 
and we can address expressions (f2l)| and (|2"T|) , as soon as we restrict the study to early times, 
so that a(t) <C a . Furthermore, using the power-law relation © for the bulk-viscosity 
coefficient, one easily finds 



Co = Co t" 



Co = z/(6ttGY 



(37) 



With the help of this expression, we can isolate two constants in the equation (|36|) . which 
finally rewrites 



5 + 



4 

3t 



X 



£2(8-1/3) 



A- 



£27-2/3 



where the constants \ and A are 



X 



t 2( S -l/3) Coq 2 

Pa a 2 



A 



2 
31 2 " 

t7 -l/3 



0, 



v s q 



(38) 



(39) 



This equation can not be analytically solved in general. Following the analysis developed in 
[24] . let us now discuss the case s = 5/6. Indeed, this case is the only of physical interest since 
it deals with the maximum effect that bulk viscosity has without dominating the dynamics, 
in view of its non-equilibrium perturbative characterization. In fact, in the collapsing limit 
as i -> oo, ifs > 5/6 the viscous term proportional to \ results to be of higher order 
and dominant. On the other hand, it can be neglected in the equation (|38|) . if s < 5/6. 
Substituting this value in the equation above, one gets the following integrable expression 



The solutions are 



<5 + 



5(t) = V 



"4 


5 


"A 2 2 " 


3 + * 


~t + 


_£2 7 -2/3 it 2 


1 X 






6 2 


C\ J n ( 


— --) + C 2 Y n 









6=0 



<■ 7 



where J n and Y n denote Bessel functions of first- and second-kind, respectively, and 

j = 7 - 4/3 . 



= -v/25 + 6x + 9x 2 /67 



(40) 



(41) 



(42) 
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These functions oscillate for t <C A 1 / 7 , while for t ^> A 1 / 7 the density-contrast solutions 
(HJ) evolve like 

5^t-Ve-x/2Ttn. (43) 

A simple analysis of the exponent of such solutions shows how it is always positive, for all 
values of the viscous parameter \, as soon as we choose the (-) sign solution. This behavior 
corresponds to a gravitational collapse, if we consider the asymptotic limit t — > oo. The 
threshold value which separates the different regimes, implying the growth of the density 
contrast, is defined by the relation t > A 1 / 7 which, using (j2"6")l , corresponds to the Jeans 
condition (TT2|) : 

v 2 s q 2 /a 2 < 6nG Po , (44) 

as in the non-dissipative case. We remark that such solutions will apply only after the 
recombination, with 4/3 < 7 ^ 5/3. In fact, in correspondence of 7 = 4/3, the solutions (T4"T]) 
show a singular behavior and the equation (|38[) requires a different treatment. 

As in the standard Jeans Model, the key value of the Jeans Mass is not affected by bulk 
viscosity, i.e., gravitational collapses for S — > 00 are addressed if 



* < Kj = \ ^ (45) 

The effect of dissipative processes is to modify the evolution of perturbations. In fact, com- 
paring expression (|43p wrt the non-dissipative behavior of growing density contrast 8 ~ i 2 / 3 , 
see (|31|) . one can show that the relation — 1/6 — x/2 — yn < 2/3 is always verified. We can con- 
clude that the effect of bulk viscosity is to damp the density contrast evolution, suppressing 
the structure formation as in the Jeans Mechanism. 



5 Concluding remarks 

The effects induced by the presence of bulk viscosity have been analyzed in two different 
cases in a perturbative scheme: the standard Jeans mechanism and the generalization in 
treating an expanding Universe background. In both approaches, viscosity has been intro- 
duced in the first-order dynamics via a power-low function of the energy density and the 
effects produced on the unperturbed dynamics have been consistently neglected in view of 
the phenomenological nature of such kind of viscosity. 

The main result, in dealing with the viscous generalization of the models proposed, has 
been to show how bulk viscosity damps the density contrast evolution maintaining unchanged 
the threshold value of the Jeans Mass. Such an effect suppresses the sub-structure formation 
in the top-down fragmentation mechanism. 

In particular, in the analysis of the dissipative Jeans Model, a new decreasing regime 
for perturbations has been found. The presence of such a behavior allowed the study of 
the top-down scheme for small and strong viscous effects. In the first case, the density- 
perturbation amplitude of a sub-structure remains substantially constant during the main 
structure collapse. On the other hand, if viscous effects are sufficiently strong, the sub- 
structure vanish in the linear perturbative regime, unfavoring the fragmentation. 

In the second part of this work, the static and uniform background solution for the un- 
perturbed evolution has been generalized by the dynamics of an homogeneous and isotropic 
Friedmann Universe. In this scheme, a Jeans-like relation has been obtained and a consider- 
able damping of the density contrast growth has been found. 
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